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Abstract 

Extending our previous results on trans-Planckian (Gs ^> K) scattering of light particles 
in quantum string-gravity we present a calculation of the corresponding S'-matrix from 
the region of large impact parameters (b ^> Gy/s > X s ) down to the regime where clas- 
sical gravitational collapse is expected to occur. By solving the semiclassical equations 
of a previously introduced effective-action approximation, we find that the perturbative 
expansion around the leading eikonal result diverges at a critical value b = b c = 0(Gy/s), 
signalling the onset of a new (black-hole related?) regime. We then discuss the main 
features of our explicitly unitary S'-matrix - and of the associated effective metric - down 
to (and in the vicinity of) b = b c , and present some ideas and results on its extension 
all the way to the 6^0 region. We find that for b < b c the physical field solutions 
are complex-valued and the S-matrix shows additional absorption, related to a new pro- 
duction mechanism. The field solutions themselves are, surprisingly, everywhere regular, 
suggesting a quantum-tunneling - rather than a singular-geometry - situation. 
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1 Introduction 



In this paper we will discuss, within a quantum string-gravity framework, the high-energy 
scattering of light particles in a variety of kinematical regimes. We will resume, for this 
task, our twenty-years old trans-Planckian S-matrix analysis [H El El E] and we will 
extend it to the situation in which, at a classical level, the initial state is doomed to 
collapse due to the appearance of a closed trapped surface [6] . 

Renewed interest in this problem stems from a growing conviction that a consistent 
quantum calculation of a collapse process leading to black-hole formation and to its subse- 
quent evaporation is the best -if not the only- way to understand the fate of the apparent 
information paradox [7J; or, better, the apparent loss of quantum coherence implied by 
black holes. Hybrid quantum arguments in terms of classical gravitational solutions do 
not meet, in our opinion, the necessary consistency requirements. It has indeed been 
suggested [8] that pure quantum states would not produce gravitational collapse even 
if the energy distribution would classically predict it to happen. And, even more dras- 
tically, it has been proposed that quantum back-reaction on the metric in apparently 
collapse-prone processes, would generate everywhere regular solutions without singulari- 
ties and event horizons [9]. Finally, topologically non-trivial (i.e collapse-like) classical 
configurations may turn out to be irrelevant in the quantum formulation of the physical 
process [10J. 

Hints of what may actually happen has to come from the actual treatment of a collaps- 
ing system in a consistent quantum theory of gravity. Unfortunately, there are not many 
candidates for such a theory. To our point of view, among them, string theory is the only 
one allowing for a treatment of the problem in perturbative, as well is non-perturbative 
regimes, despite the fact (or perhaps, because of the fact) that it is not, to start with, 
a general relativistic theory describing space-time dynamics. Strings can only be consis- 
tently quantized in appropriate backgrounds, those that do not introduce two-dimensional 
Weyl anomalies. As in our previous papers, we will study the scattering process in D = 10 
superstring theory in Minkowski space-time -after compactifying n dimensions on string- 
size tori- and look at possible (perhaps even approximate) interpretations of the results 
in terms of an effective metric best describing the quantum process. 

Let us briefly recall our (ACV hereafter) approach and results. Scattering of two 
massless strings (e.g. of two gravitons) was considered at centre of mass energy 2E = 
t/s ^> Mpi anc k and impact parameter b in (^-dimensional Minkowski spacetime, where 
d = 10 — n. In this paper we shall work in d = 4, but we expect that the extension to 
d > 4 will not present major problems. We shall also focus on a regime in which string-size 
effects are relatively small, while the gravitational interactions can be strong. In order to 
define this more precisely, let us recall that, in string-gravity, the fundamental scale is the 
string length X s = y/a'h, in terms of which the Planck length and the Newton constant 
are expressed as Xp = \J Gh = gX s , where g <C 1 is the string loop expansion parameter, 
assumed to be small. 

On the other hand, at very high energies y/s, the gravitational (Schwarzschild) radius 
R = AGE = 2G\fs plays an important role. In our small-coupling, high-energy regime, 
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defined by Gs » H, R is much larger than Xp, but can be smaller or larger than 
X s , because the ratio R/X s = g^jGsjh involves the small coupling constant g. As a 
consequence, there are three distinct regimes according to which one of the three length 
scales b, R, or A s exceeds the other two. 

If b ^> R, X s one deals with small deflection-angle scattering. This is well decribed by a 
leading eikonal approximation with small string-size and classical corrections correspond- 
ing to the expansion parameters (X s /b) 2 and (R/b) 2 , respectively. The former are quite 
easily taken into account pQ and can best be interpreted [11] as string excitations due 
to the tidal forces induced on each string by the effective (Aichelburgh-Sexl) shock-wave 
metric produced by the other string. They have been analyzed by ourselves in the past 
[T] and, most recently, in [T2] 

In the regime A s 3> b, R, also investigated through fixed-angle scattering [T3] , string 
effects soften gravity according to the generalized uncertainty relation [TH [2] 

h 

Ax > — + a'Ap > X s . (1.1) 
Ap 

As a consequence, the minimal observable size of the system is A s itself, which exceeds 
R, and classical gravitational collapse conditions are never met. It is possible, however, 
to push the analysis of this regime towards its boundary X s — > R > b, which should 
correspond to the threshold for black hole formation threshold ~ M s g~ 2 ~ M P g~ x . One 
finds JT5J, [1] that, even if no black-hole is formed, the final state, in the energy region 
Mp < E < -^threshold starts to vaguely resemble that of an evaporating black hole of mass 
y/s with typical final momenta of order Mp/^/s ~ h/R. In other words, a precocious 
black-hole-like behaviour is found to occur even below the expected threshold for their 
actual production. 

By contrast, for R > X s , new semiclassical phenomena take place. They extend beyond 
the impact parameter at which string fluctuations, including those due to diffractive 
excitations, are large. This is the regime that we attempt to treat in this paper, for 
various values of the impact parameter b > X s of the colliding strings. The interesting 
region is the one in which b approaches R from above and possibly goes below it, a 
situation in which, classically, a gravitational collapse would take place [5]. A general 
framework for describing this most difficult regime was proposed in [16], where the S- 
matrix was connected to properties of a classical solution at past and future null infinity 
(the so-called Bondi masses). However, in spite of its conceptual appeal, going beyond 
the leading eikonal in that formalism has proven prohibitively difficult. Here we shall 
use instead our key observation pQ that, because of the softness of multi-loop string 
amplitudes, the S-matrix exponentiates in terms of an eikonal function of order Gs/h 
which, in turn, can be expanded in powers of R 2 /b 2 for b > R. The outcome has a 
diagrammatic interpretation that can be encoded into an effective action. If string effects 
are neglected, that action agrees with Lipatov's effective action [TF] (see also [18]) and 
reproduces [4] the previously computed leading-order correction to the eikonal [3]. 

Therefore, the effective Lagrangian that we investigate in this paper is motivated by 
our string-gravity expansion, even if it does not contain explicit string corrections. It 
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is a function of appropriate components of the metric h^ v {x) which are apt to describe 
the high-energy regime, interact among themselves via the effective coupling R/b, and 
are coupled to sources provided by the energetic scattering particles. The solutions of 
the (nonlinear) lagrangian equations provide an effective metric (which appears as the 
outcome of quantum backreaction effects) in terms of which the action, and thus the S- 
matrix, is expressed and computed. The unitarity of the approach implies an absorption 
in the elastic amplitude due to particle production, whose inclusive properties (spectra, 
correlations, etc..) may be analyzed. 

Within this framework, we treat here the region b > R and we attempt to tackle the 
most interesting region b < R where still one should be able to compute the scattering 
amplitudes and the effective metric. We cannot claim to have fully achieved that goal, 
but we have progressed pretty far towards it, even beyond expectations. We thus believe 
that the analysis of our results should provide at least some hints as to whether - in this 
consistent quantum approach - there is any sign of a trapped region or event horizon and 
what is the "unitary evaporation" that is produced without loss of quantum coherence. 

The rest of the paper is organized as follows: 

In Section 2 we recall the eikonal expansion, the form of the first correction to its 
leading term, and the effective action that should generate through its tree diagrams 
the higher order corrections. We then define a somewhat simpler problem in which a 
"rescattering" term, the related "double-diffractive" string excitation, as well as one of 
the emitted-graviton polarizations, are neglected. The above approximations will be used 
throughout the paper, although, in Section 6, we will give an educated guess on how the 
second (infrared-sensitive) polarization can be included. 

In Section 3 we discuss the axisymmetric case in which the field equations become 
ordinary differential equations. We first consider a class of analytical solutions to the field 
equations in the case of point-particle collisions at b = 0. We are aware of the fact that 
this is a most difficult regime for justifying some of our approximations, in particular the 
neglect of string corrections, which should be restored later on. Nevertheless, this class 
of b = solutions - which is surprisingly simple and robust, but quite non-perturbative 
- turns out to tune up the discussion on the boundary conditions to be set in order to 
match perturbation theory at larger values of b > R. Furthermore, they are likely to be 
essential for the overall interpretation of the problem for A s < b <C R. As an amusing 
digression we will also discuss here the central collision of two extended sources (taken to 
be two identical homogeneous disks of radius S for simplicity) where the problem can be 
solved analytically and shows the existence of a critical ratio S/i?. This case could also 
be a way to represent string-string collisions with E ~ \ s . 

In Section 4 we turn to the case of generic values of b where, within some technical 
approximations, we are still able to solve the problem analytically. We find that, while at 
b~^> R the perturbative expansion is qualitatively correct, the expansion diverges at some 
calculable critical value of b = b c ~ R. We also discuss possible ways to define solutions 
below b = b c . 

In Section 5 we reformulate the problem in momentum space as a set of integral 
equations lending themselves to an iterative solution. We reach conclusions that are in 
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very good agreement with those obtained in the position space approach of the previous 
section: in particular, the iterative solution only converges above a critical value of b/R. 

In Section 6 we turn our attention to the construction of a unitary S'-matrix, to the 
properties of the final state and to the expectation value of the metric in that state. Once 
more, the properties of the final state appear to resemble those due to an evaporating 
black hole as we approach b c . 

In Section 7 we describe (and try to interpret) our proposal to define the scattering 
amplitude and effective metric for b < b c , which is based on the analysis of properly iden- 
tified complex- valued solutions of the field equations. Finally, in Section 8, we summarize 
our main results and give a brief outlook. 



2 Eikonal Expansion and effective action 

In string-gravity, ACV found that the S-matrix in the impact parameter representation 
has an eikonal form where the eikonal operator can be expanded in the parameter R 2 /b 2 . 
For b » R > X s , the eikonal resums all powers of Gs due to multigraviton exchanges, as 
follows (see also [13]): 

Os L 

S(b,s) = exp2i6 (b,s), 5 = — log(-), (2.1) 

where L is an infrared cutoff related to the well known infinite Coulomb phase. String 
effects in this region are simply taken into account [lj by an operator shift of the impact 
parameter variable 

8 ^ 8 = 8 (b + X u -X d ,s), (2.2) 

(where an average is performed over the closed string position operators X u ,Xd) and 
give rise to the diffractive string excitation and fluctuations mentioned before. The soft 
behavior of multi-loop string amplitudes is itself responsible for the dominance of eikonal 
iteration in the results (12.11) . (12.21) . In fact, Eq. fj2.lH can be interpreted (fig. (pQ)) as a 
multiple scattering series, in which the (small) deflection angle 9 = 2R/b - corresponding 
to a possibly large momentum transfer t ~ Gsh/b 2 - is built up by many graviton-exchange 
processes of small momentum transfer, of order t s ~ (H/b) 2 <C (h/X s ) 2 . 

When R/b becomes sizeable, ACV found that the eikonal can be expanded in a power 
series in R 2 /b 2 and, possibly, X 2 /b 2 . The terms (R 2 /b 2 ) n of such a series are in correspon- 
dence with connected tree diagrams interacting with the colliding strings via the exchange 
of 2n (reggeized) gravitons, as shown in fig. (j2J). Besides the one-loop correction 8\(b, s), 
ACV found that the lowest term in such a series is the so-called H-diagram of fig. ([3D, 
contributing at two loop level to the real part of 82 

This extra contribution to the phaseshift modifies the Einstein deflection angle of energetic 
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Figure 1: The leading eikonal rescattering series. Crosses denote on-shell propagation 



(massless) particles in the form: 



6 r i R R 2 



(2.4) 



Further terms in the expansion are expected, and will be calculated here in a framework 
to be defined shortly. 
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Figure 2: Diagrammatic structure of 2n-loop irreducible contributions to the eikonal 

ACV found also inelastic effects contributing to the imaginary part of the phaseshift. 
There is an infrared divergent contribution, connected with soft graviton bremsstrahlung, 
not explicitly discussed here, and a finite part Im5 2 = 2 log s Re^/V, connected with hard 
graviton emission, which will be generalized in the following. 

2.1 H-diagram: Amplitude and emission field 

In order to extract from the H-diagram in fig. [3] a computational method for higher 
orders in R 2 /b 2 of the eikonal operator, we have introduced in ref. [I] an effective action 
approach. An essential ingredient in it are the high-energy graviton [20] and string [21] 
emission vertices, which lead to the emission amplitude that we now recall. 

While the leading eikonal exponential is generated by the exchange among the col- 
liding particles' sources of an arbitrary number of (longitudinal) gravitons, the ~ R 2 /b 2 
correction is represented by a graviton being emitted by two exchanged ones and then 
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absorbed by two others. We shall see that this intermediate graviton, in its transverse 
polarization, will play an important role in the scattering process. To this purpose, let 
us define the following two independent transverse-traceless polarization tensors for a 
graviton of momentum k ( bold-case notation referring to transverse momenta): 



uv 
-TT 



(J, v 



L C LJ > C LT 



fU/ 3 



where: 



e1 + ie^ = ^(k) = (— ,ie,— ), 



25) (i,j = TT, LT) , (2.5) 

(2.6) 



and e is the unit polarization vector transverse to k. 

At high energies, on the basis of the vertices in [20], [21], the graviton emission am- 
plitude of fig. ([3]) takes the form {h = 1) 



K 3 S 2 
1 2 

2/cV 



Re[(k 2 kf - \k l \'\k 2 \ I )r 1 ' M r 1 



(2.7) 



k 2 



sin Q\ 2 e^ T — sin B\ 2 cos 0i 2 e^ T ] 



where we have defined k 2 = 8nG and k = ki + k 2 = k 3 + fc 4 . From eq. (12.71) we derive 
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Figure 3: Regge-Gribov H-diagram, yielding the first subleading correction to the eikonal; 
dashed (wavy) lines denote exchanged (emitted) gravitons corresponding to the fields a 
and a (h Q or <p Q ) 



the imaginary part of the H-diagram [3] 
lmA H (s,q 2 ) _ Y 

S 167TS 2 



d[k 1 ]d[k 2 ]A^(l,2)A*(3,A), 



(2.8) 



where d[k] = d 2 k/(2ir) 2 , Y = logs and q = k 2 — k 3 . We then obtain the real part by 
a dispersion relation which amounts to multiplication by ir/2Y, and finally the impact 
parameter amplitude by a Fourier transform □: 

(SttGQV f e-«« . 2 . 2 sin 29 12 sin 2fl 34 , 
Re 5H[b,s) = — / d[fci]d[fc 2 ]%]-^(sm 9 12 sm 6 34 H - ) . (2.9) 



1 We use conventions in which AsS(b, s) = J d[q] e zbq A(q, s) and ImA(0, s) = soto 
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By then introducing a fc 4 -mtegration through a delta-function 5(ki + k 2 — k 3 — fc 4 ) we 
can rewrite eq. (12.91) in two different convenient forms. The first, to be used in sees. (5) 
and (6), uses the simple identity 2q = {k 2 _ &i) + (&4 — ^3) to yield the factorized form: 

Re6 H (b,s) = ^lm6 H (b,s) = ^Gs(nR) 2 J d[k] (\h TT (k)\ 2 + \h LT (k)\ 2 ) , (2.10) 

where we have introduced the contributions of the TT and LT polarizations in eq. (12.51) : 

h TT (b, k) = 16tt 2 J d[kl j^ k2] 5(k -k x - k 2 ) exp (ibk 2 )) sin 2 9 12 , (2.11) 

h LT (b, k) = 16vr 2 J < ^^^-5(k-k 1 -k 2 )exp(ibk 2 ))sme 12 cose 1 2. (2.12) 
Alternatively, we can rewrite eq. (12.91) as an as-space integral 

Re 8 H {b, s) = ttGs^^- J d 2 x(\h TT (b, x)\ 2 + \h LT (b, x)\ 2 ) , (2.13) 
where the analogous contributions of the TT and LT polarizations in positions space read: 

J (k 1 + k 2 y n 2 \b-x\ 2 

h LT (b, x) — 4 / TTr^—rrb sin cos ex P i(k 2 b - kx) . 
J (Ki + k 2 ) 

Here we note that Iitt in eq. (12.141) has a simple expression, which curiously reproduces 
the form of its Fourier transform and has a 1/cc 2 behaviour at large distances. On the 
other hand, Jilt is more involved and shows a l/\x\ = 1/r behaviour for r ^> b that, in 
turn, produces the well-known logarithmic infrared divergence in eq (12.91) . due to graviton 
bremsstrahlung. Indeed, one can show [3] that the LT polarization is responsible for the 
Weinberg current [22]. The corresponding infrared behaviour was discussed in detail 
in ref. [3], where a subtraction in dimensional regularization was performed in order to 
obtain the finite result in eq. (12.31) for 5 2 . On the other hand, here we are interested in the 
possibly collapsing energy, not in the one which is peripherally radiated away. Therefore, 
in most of the following, we will subtract the LT polarization altogether, by restricting 
ourselves to the TT one, which is IR safe. 

Let us note some interesting properties of the /i-fields just introduced. By defining 
z = x 1 +ix 2 and d = d/dz, we find that the complex combination ho = IiTT + ihiT, by the 
i sin 6 12 exp(— i6i 2 ) form of the couplings in eq. (I2.14p . satisfies the differential equation 

2\d\ 2 h (b,x) = 4(d 2 a d* 2 a -\d\ 2 a \d\ 2 d ) = 1-_J_ , (2.15) 

7T Z* z (b — Zf 

where we have defined the fields a and a - to be related to longitudinal gravitons - by 

1 \z\ 2 _ 1 
a (2;) = -— log(— ), a = a (b-z); \d\ 2 a = --5(x) . (2.16) 
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One may also double check that the expression for h^r given in (12.141) satisfies the real 
part of (I2TT51) i.e. 

V 2 Re/* = Wt = * i 2 ^'^ - l) . (2.17) 

n 2 x 2 (b — x) 2 \ x 2 (b — x) 2 I 

Furthermore, one can define the Fourier transform of eq. (12.71) : 

= 2kV J (sin 2 la e^ T - ~ sin 26 12 e%\ expt(k 2 b - kx) , (2.18) 

and an effective gravitational field related to the H-diagram. This can be written in terms 
of the h field as 

^- = ^(x) = ^Re[^h (b,x)}, h%=??ho(b,x) = ^ V y2 ° l °V (2.19) 

where we have promoted the polarizations e l to operators in a?-space. The result is better 
rewritten by introducing a (generally) complex scalar field <fi such that @. 

h = 4|<9| 2 = V 2 0; hij = (SijV 2 - didj)(f> . (2.20) 

In this language, restricting to the IR safe polarization means considering Re0 only, or 
the field to be real. By replacing in eq. (12.131) the corresponding expression of hxr m 
eq. (12.141) . we obtain: 

[Re5 H (b,s)] TT = kGs { ^- J d 2 z(V 2 Re0 o ) 2 = Gs^ , (2.21) 

to be compared to the full result Re£ 2 (&, s) = Gsf^ of eq. (12.31) . the difference being due 
to the (neglected) LT polarization. 



2.2 The reduced effective action 

Further terms in the i? 2 /6 2 -expansion are obtained by considering both multi-H diagrams 
combining multiple emissions (fig. (jl|)) and rescattering diagrams in which emitted gravi- 
tons reinteract by a longitudinal exchange (fig. (E])). Here we shall limit ourselves to the 
first class of diagrams, that will be treated to all orders, while the second class - which 
starts at order i? 4 /6 4 - will be briefly discussed in the next subsection. 

Multi-H diagrams are described by a (reduced) two-dimensional action introduced in 
ref. jl] in which the longitudinal fields a and a and the mostly transverse field h occur, 
by generalizing the contributions of eqs. (I2.15P and (I2.16p . By restricting ourselves to the 

2 Here the field 4> has a more convenient normalization, for writing the action, than that of ref. [4]. 
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Figure 4: Regge-Gribov double- H diagram, contributing at order R 4 /b 4 to the eikonal; 
dashed (wavy) lines denote exchanged (emitted) gravitons 



IR safe polarization, h and are now supposed to be real valued and the reduced action 
takes the form 



2ttGs 



a(6) + o(0) - i y d 2 a;VaVa + ^ d 2 z:(-(V 2 0) 2 + 2ftV 2 



-V 2 ft = V 2 a V 2 a - ViVj-a V^-a . (2.22) 

Here the longitudinal fields a(x), a(x) interact with point-like sources placed at x = b 
and x — 0, respectively. The field <fi is generated by the current Tl, which is defined 
by a generalization of eq. ( 12.151) and expressed through vector derivatives, as it's more 
appropriate for real-valued fields. Note in ( 12.22ft the appearance of the effective coupling 
R 2 = 4G 2 s, which controls the dependence of the solutions on the expansion parameter 
R 2 /b 2 . 

The coupled lagrangian equations derived from eq. (12.221) read 

V 2 a + 25(x) = 2(vr J R) 2 (V 2 a V 2 0- ViV^a ViVj-0), a(x) = a(b - x) 

\7 2 H = V 4 = - (V 2 a V 2 a - ViV 3 -a ViV^a) . (2.23) 

It is soon apparent that TC — h — V 2 </> on the equations of motion, so that the "on-shell" 
action reads 



72 i\2 



4%^ = a(6 2) + a(o) _ 1 f cfxVWa + ^f [ d 2 x{V\, 
2nGs 2 J 2 J 

= 2a(b)+I a + I (t> . (2.24) 

By performing the first perturbative iteration, we find thet a, a and h reduce to the 
expressions in eqs. (12. 16[) and (I2.15p . while the low order contributions to the action 
become 

/(°) = -oo(6), J« = -2ax(6) = -Uf , (2.25) 
so that, by collecting all terms, they partly cancel and finally yield 

A(b,s) = 2nGs(a (b) + J< 1} ) = 2nGs a (b) + 2Re a H = 2Gs(log \ + , (2.26) 
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thus reproducing the perturbative result for the TT part. 

The framework defined by eqs. (12.221) and (12.231) is the one we shall analyze in detail 
in the following, including its features in the gravitational collapse region in which the 
effective coupling R may exceed the physical size of the system b, larger than the string 
size X s . Besides lacking the LT polarization - that we have argued to be related mostly 
to an infrared, peripheral phenomenon - the above model lacks rescattering and string 
effects, whose form we shall briefly recall in the framework of an effective action [T71 [TS] 
which includes the dependence on the light-cone variables x ± = x° ± x 3 . This will allow 
us to relate the reduced action to a shock-wave solution having a particular form in the 
x + , x~ plane. 

2.3 Shock-wave interpretation 

The reduced action model introduced above can be viewed as a particular limit of Lipatov's 
effective action [TTJ [18], which in turn provides a formal description of the diagrammatic 
series in Fig. (T5]) in the limit in which all subenergies among emitted gravitons are large 
and string excitations are neglected. The corresponding Lagrangian contains, besides a 
field $ (related to the previously introduced </>), the longitudinal fields h ++ and h which 
are similarly related to a and a and are coupled to the external sources of the impinging 
particles (gravitons or strings 

In the effective action framework, the elastic S-matrix of the tree diagrams in fig. ([2]) 
is given in terms of the classical solutions of the lagrangian equations of motion as 

S(b,s) = ex P [^(/C)]; (2-27) 
A(h ++ , h—,$) = J d 4 x(C + C e + C r + T ++ h ++ + T-Jk~) , 

where &(x + ,x~ ,x) generalizes 4>(x) to four dimensions, and 

T__ = kE8{x~)8{x), T++ = kE5(x + )5(x - b) (2.28) 

represent (up to an unconventional but convenient factor of 2k), the energy-momentum 
tensor of the colliding particles. The lagrangian consists of a kinetic term 

C = -d*~h ++ dh-~ + 4<9+<9* 2 $<9_<9 2 $* , (2.29) 

where the longitudinal fields have a mostly transverse propagator and the (complex) $ 
field a mostly longitudinal one, of a graviton emission term 

C e = K(J\d\ 2 ^* + J*\d\ 2 <S>); \d\ 2 J = [d* 2 h ++ d 2 h— - \d\ 2 h ++ \d\ 2 h— } (2.30) 

related to the reduced one in eq. (12.221) . and, finally, of a rescattering term 

C r = K(h ++ d* 2 &d + 2 d 2 <& + h~-d 2 <&*dJd* 2 §) , (2.31) 

3 As better explained in sec. 6, our fields are related to the usual metric components by: h^dx^dx" = 
ds 2 - rj^dx^dx" = 2 K (h ++ (dx+) 2 + h—(dx-) 2 ) + (k/4)(c£J ARe<I> - ej^AIm^)dx^dx u 
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which is supposed to take into account the rescattering diagrams of fig. This term is 
quadratic in $, and is likely to play a role when the latter is large. 



Figure 5: Typical Regge-Gribov rescattering diagram, contributing at higher orders to 
the eikonal; dashed (wavy) lines denote exchanged (emitted) gravitons which have mostly 
transverse (longitudinal) propagators. 

Here we do not treat the action (I2.27P in detail, but we would like to discuss a couple of 
important points. First of all, it was shown in j3] that the reduced action and equations of 
sec. (12. 2p correspond to a shock-wave solution of the present lagrangian equations without 
rescattering terms, of the form 

h ++ = ftv^s 6(x~)a(x), h = 5(x + )a(x); (2.32) 

/Tj S 

$ = — <d(x + x~)(f)(x) , 

where now the longitudinal fields a, a and the transverse field <fi appear as profile functions 
in front of the x + , x~ dependence. Note that, while the longitudinal part is of Aichelburg- 
Sexl type, the transverse part has support inside the whole light-cone. This propagation, 
of retarded plus advanced type, corresponds to the principal value part of the Feynman 
propagator and is appropriate for the real part of the amplitude B 

The shock-wave interpretation of the reduced action framework allows to embed it in 
spacetime, and to provide, in particular, the effective metric produced by the solutions of 
the lagrangian equations. The explicit x ± -dependence in eq. (I2.32p allows to calculate the 
longitudinal components of the metric induced by the $ field by generalizing eqs. (12.181) 
and (I2.19P to x ± -space and by using the longitudinal components of the TT polarization 

TT _ ®+ TT _ 8- TT _ 1 / oo\ 

Taking into account the different normalization of the metric components mentioned be- 
fore, this procedure leads to the following expression for the induced metric in terms of 



4 Restoring the transverse propagation corrections to the <!> field solution amounts to effectively cutoff 
the wave inside the light-cone around x + x~ < R 2 without modifying the wave-front, and thus the 
derivation of the reduced action. 
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a, a and a real- valued 0: 



ds 2 



dx + dx-{\ - -(vr J R) 2 6(x + x~)V 2 0) + 2vri? (a(^)5(x")(rfx~) 2 + a(z)<5(x + )(cfe + ) 2 ) 
(tt j R) 2 V 2 (|x+|5(x-)(c/x-) 2 + |x-|5(x+)(rfx+) 2 ) + ds 2 T 



1 



4 




|dz| 2 + (7TR) 2 Q(x + x-)(2\d\ 2 (f) \dz\ 2 - <9 2 dz 2 - <9* 2 dz* 2 ) 
\dz\ 2 + (ttR) 2 Q(x + x-) (5 ij V 2 - ViVj) dx l dx j . 



(2.34) 



It is easy to check that the perturbation of the metric proportional to is transverse and 
traceless: it has exactly the form of the TT polarized gravitational field introduced in sec. 
(2.1), and is meant to describe the intermediate h field contributing to the real part of 
the amplitude. A discussion of the features of the above effective metric is postponed to 
sec.fllj), when explicit solutions will be available. 

Note finally that we shall not consider in the following rescattering and string contri- 
butions to higher orders of the eikonal expansion. This is an acceptable approximation 
for the small regime R < b, but is likely to be insufficient when looking at distances 
smaller than R and approaching the string length, when becomes large (cf. sec. (J3J)). 
Note that in this region rescattering and string effects are probably intertwined, because 
of the eikonal couplings occurring in the rescattering vertex of eq. (12.311) . They produce 
factors of k + k~ in the numerator of the corresponding diagrams and thus make the lat- 
ter formally divergent, by emphasizing the role of large intermediate masses. Therefore, 
string excitations are required in order to regularize the sum over intermediate masses, 
and become non negligible. We have here a situation similar to normal diffractive exci- 
tation of string massive states by initial particles, where string corrections are taken into 
account, eventually, by the simple shift in eq. (12.21) . In the present case it is the interme- 
diate graviton to be excited in a sort of double-diffractive string excitation. We similarly 
hope that the present effects will turn out to be calculable, perhaps by introducing, in 
the framework of sec. (12.21) . the Regge-graviton string emission vertices of ref [21]. 

3 Axisymmetric solutions 

3.1 Particle-particle scattering at b = 

We start considering the reduced action of sec. (12.21) in the complementary region to the 
perturbative one we started with in sec. (12.11) . by taking X s <^ b <^ R. Since b > X s we shall 
not consider string corrections explicitly, even if the string has played an important role 
in assessing the validity of the model. This means - since R is the only explicit coupling 
being considered- that we actually take the 6 = limit of a head-on collision. This is a 
fully non-perturbative regime whose interpretation requires a non-trivial matching with 
the perturbative regime that we shall discuss in the following sections. 

In the 6 = limit, we can look for axisymmetric solutions a = a and which are 
functions of r 2 = x 2 only. Surprisingly, the nonlinear equations (12.231) take a particularly 
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simple form that will allow a complete treatment of their solutions. Indeed, by setting 
a = da/dr 2 etc.., they read, for r ^ 0, 



_d_ 

dr 2 



[r 2 d(l - (2vr J R) 2 0)] = 0, (3.1) 
d r „ 1 d 



nr y^<f>)] + -—{r 2 a 2 ) = 0, (3.2) 

thus providing, by inspection, two constants of motion with respect to the "time" variable 
r 2 . The first one in eq. ( 13. ID is fixed by the Gauss theorem on the delta function in (12.231) 
to be — 1/27T and the other vanishes by the same token, so that we have 

r 2 d(l - (2vr J R) 2 0) = C X = ~, r 2 (r 2 + \a 2 ) = C 2 = . (3.3) 

Z7T 2 

It is now convenient to introduce the function: 

p(r 2 )=r 2 (l-(2nRf<j>), (3.4) 

which has dimension of a squared length, and embodies the effect of the (transverse) 
emission field which will play an important role in the following. By expressing <fi and a 
in eq. (I3.3P in terms of p we have 

a ^ = ~ (3 - 5) 

Kr 2 )-$- 2 =0, p 2 + — = C s , (3.6) 
2p 2 p 

where the constant C3 will be determined by requiring consistency with the perturbative 
expansion for r 2 3> R 2 . Indeed, in this limit, the system of eqs. (12. 23ft reduces to eqs. (12.161) 
and ( 12. 151) thus showing, by (12.141) . the large-r behaviour 

1-1 R 2 

a(r 2 ) = a(r 2 ) ~ logr 2 ; r 2 <p ~ — -logr 2 , p(r 2 ) ~ r 2 log r 2 , (3.7) 

2tt Six ' 2 

which implies p — > 1 at large distances and thus C3 = 1 in eq. (13.61) . yielding finally 

P 2 + - = 1, a = (3.8) 
p lixp 

We conclude that our nonlinear problem reduces to the classical Coulomb problem for 
the "radius" p{r 2 ) at "time" r 2 . Since the "Coulomb potential" in eq. (13. 8p is repulsive, 
we can say from start that, coming from a positive p(r 2 ) ~ r 2 at large distances, the 
generalized radius p(r 2 ) > R 2 will never vanish during r 2 -evolution, even at r 2 = 0. This 
means that d in eq (13.51) is not singular at r 2 = - contrary to its perturbative behaviour 
- and that, by eq. (13. 3j) . (27ri?) 2 ~ — p(0)/r 2 must be singular instead, a feature which 
is non-perturbative as well. 
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The above non-perturbative behaviour of a and <fi is somewhat puzzling. It means 
that in the first eq. (13. 3p the outgoing flux of Va is traded for that of V0, with increase 
of the <p field at small distances and a non-vanishing value of r 2 0, due to p(0) 7^ 0. It 
is caused by the large curvature p for small values of p which forbids p(0) = for non- 
negative, real-valued solutions. On the other hand, the condition p(0) = is not only a 
property of the perturbative behaviour, but appears to be required in order to avoid or 
to treat properly a boundary term at r = in the reduced action, as discussed in the 
Appendix. We shall see that only for b > R we will be able, in the next section, to meet 
that condition for real-valued solutions, thus obtaining a nonsingular 0. Alternatively, we 
can give up the reality requirement suggested by the metric interpretation, and look for 
complex 6 = solutions. We shall motivate and explore such possibility in sec. 7. 

The explicit solutions of eqs. (13. 8p for p and a = a are obtained by standard methods 
in terms of a hyperbolic angle x( r2 ) as follows 

p(r 2 ) = i? 2 cosh 2 x(r 2 ); r 2 = R 2 (x + coshx sinh x — Xo ~ cosh xo sinhxo) (3.9) 

a{r 2 ) = f f 2 = ~{X{L 2 ) - X{r 2 )) , (3.10) 

2n J r 2 p(H) 7T 

where xo — x(0) is the arbitrary value of x at the origin. Its presence is not surprising, 
because we have set only one boundary condition, providing matching to perturbation 
theory at large distances (p(°°) = !)• The explicit solution (13. 9p shows that the boundary 
condition p(0) = can only be met by a complex value of Xo (like Xo = —in/ 2) - a case 
that will be discussed in detail in sec. 7. For the moment we treat Xo as a f ree parameter, 
even if a way to determine a real value for it will be discussed in the next section, as an 
alternative to p(0) = for b < R. 

Since xi 7 " 2 ) is a monotonically increasing function, there are two kinds of real- valued 
solutions, depending on the sign of Xo- If Xo > 0, p{r 2 ) increases monotonically to 
p ~ r 2 at large distances, while for xo < p decreases first to its minimum p = R 2 - 
corresponding to x = - and then increases to 00. The scale of the large-r 2 behaviour 
of p is itself dependent on xo- Indeed, a simple iterative evaluation of eq. (13.91) yields the 
more detailed behaviour 

R 2 4r 2 

p ( r 2) ^ r 2_ log _^^. ^2^ = #2 ( 1 + 2 +sinh2 ) ( 3n) 

2 r 2 (xo) 

1 4r 2 1 4r 2 

log-^; ^ -— log 2 ^--; (3.12) 



37r 2 r 2 ?" 2 (xo)' 167T 2 ^ 2 (xo 

If L 2 R 2 , 4r 2 
27 V° g ^ + Yr 



*(r 2 ) = -(^g^ + —\og—), (3.13) 



which is actually valid for any value of Xo- 

Let us remark that the arbitrary constant occurring in the integration of (13. 8p for a(r 2 ) 
has been traded, in eq. (I3.10p . for a scale L, defined by a(L 2 ) = 0, that will play the role 
of infrared scale, as in eq. (12.161) . The latter is then fixed to be the same IR cutoff L ^> R 
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needed for the evaluation of the action below, so that a(0) carries the large logarithm 
X (L 2 ) ~ logf . 

The reduced action of eq. (I2.24p can be explicitly evaluated on the axisymmetric 
solutions, and takes the form 

M[) - S) - a(0)+a(0)-2vr t dr 2 r 2 a 2 + 2ir(27rR) 2 [ dr 2 (3-14) 



2itGs v ' v ' J J dr 2 

= 2a(0) + I a + I <t>) 

where the integral I a is IR divergent, but can be made finite by combining it with a(0) as 
follows 

1 F L2 dr 2 r 2 1 

h + a(0) = - / —(!--) = -(!- e~ 2 n, (3-15) 
2tt J p p 2tt 



dr 2 (l - p) 



exp(-2xo) 



2-nR 2 J v r/ 2tt 

where the finite L — ► oo limit of J a + a(0) and have been easily evaluated by use of 
eqs. (13.91) and (I3.10p . After simple algebra we obtain 

1 2L 1 

.4(0, s) = Gs (2vra(0) + 1) = 2Gs ( X (L 2 ) - Xo + -) ~ 2Gs (log — - Xo + -) , (3.16) 

which shows a simple additive dependence on xo and the expected IR divergent Coulomb 
phase. The latter is unrenormalized, and can be factorized away in the S-matrix, as usual. 

Inserting the solution (13.91 13. lOj) in Eq. (I2.34j) we obtain an explicit expression for 
the metric where we note the appearance, besides that of A<f>, of the field <fi that - if 
p(0) 7^ - generates the behaviour r~ 2 in some metric components. Further comments 
on this issue and a discussion of the effective metric will be given after having extended 
our analysis to generic values of b. 

3.2 Central collision of two homogeneous beams 

An interesting case that can also be solved analytically but, unlike the previous one, 
contains a tunable parameter is that of the central collision of two homogeneous, finite- 
size beams of massless particles. The point here is that our effective-action method should 
retain its valididity even when the point-like sources are replaced by smooth (null) energy 
distributions on the transverse plane. A particularly simple case is that of two circular 
homogeneous beams of radius £ (area vrS 2 ), each one carrying a total amount E of 
energy, and undergoing a head-on collision. The problem is again axisymmetric and is 
characterized by the dimensionless parameter: 

^ = 4— = 4vrG e S = 4, (3-17) 
E S 2f ' K J 
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where e is the energy density per unit area and / = (87rGe) _1 is the focal distance for 
geodesies impinging on the beam-shaped shock wave (see e.g. [23]). 

At classical GR level, the problem of determining when a closed trapped surface (CTS) 
is produced by the collision led to the conclusion [6] that a CTS forms when the above 
ratio exceeds a critical value, for which an upper limit was established: 

777 <!■ ( 3 - 18 ) 



W CTS CTS 

It is quite interesting to investigate this problem within our present quantum approach. 
It is straightforward to adjust our effective action equations (13.31) to this new situation. 
They get simply modified as follows: 

r 2 d(l - (27ri?) 2 0) = -i-0(r - S) - ^0(E - r), r 2 + ^a 2 = , (3.19) 
from which, using again p as defined in (13. 4ft . we obtain 

p(0 = |^(r - E) + J^(£ - r) . (3.20) 

In other words, the equation for p is unchanged at large r > £ but is strongly modified 
(though in a continuous way) for r < S. At small r the equation has a regular solution 
with p(0) = and a nice, analytic expansion around r 2 = 0, which can be computed after 
inserting some value for p(0). This solution, however, should match the one from r > £ 
at the r = £ boundary. In this latter solution one has, as before, 




p = +\l =S> p(S 2 ) = +W1 - -7=5r = tanh Xs , (%s = x(£ 2 )) , (3.21) 



and thus the initial condition on p(0) has to be chosen so as to satisfy ( 13.211) . This, 
however, turns out to be impossible if: 

S 2 p(S 2 ) = S 2 tanh xs < p(S 2 ) = R 2 cosh 2 xe , (3.22) 

simply because the concavity of the p-curve (due to p(r 2 ) > 0) will prevent such a curve 
to pass through the origin, similarly to the previously discussed 6 = case. Such a simple 
concavity argument gives an upper limit on the value of R/H for which the condition 
p(0) = can be imposed for real-valued solutions. Equivalently, it gives a lower bound 
for the critical value of S, E c , for that to happen. It is easily computed to be given by 
{R/E) c < 2 1 / 2 3- 3 / 4 ~ 0.62. 

On the other hand, we can also provide an upper bound on E c , and thus prove the 
existence of two distinct regimes, by noting that, by eq. (13.201) . 

< ^ - = Jl [ *^ < ^-^ ry (3.23) 
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By replacing in (13.231) the values of the solution for r > £, we get the relation 



< E 2 tanh xs - R 2 cosh 2 xs < 



R 2 



(3.24) 



4 tanh 2 xt, ' 



where we realize that the r.h.s. is of order R 2 , without any particular enhancement when 
E > R increases, so that a solution for RfE c can be found. Therefore, for E > E c , the 
condition p(0) = can be met by real-valued field solutions. 

In order to get a more precise estimate of (i?/E) c we have solved numerically the 
differential equation ( 13.201) and looked for a critical value above which it is no longer 
possible to impose the condition p(0) = 0. The result of such an analysis gives: 



in line with the classical CTS-bound ( 13.181) . 

The above discussion suggests a (loose?) correspondence between "untrapped" clas- 
sical GR solutions and our real-valued field solutions that satisfy p(0) = and match 
perturbation theory at large distances. In either case such solutions cease to be available 
when the beam size is smaller than some critical radius below which gravitational trap- 
ping occurs on the classical side, and the small-r 0- singularity develops on the other. We 
shall find a similar phenomenon in the case of particle-particle scattering at b > 0, to be 
discussed next. 

4 Extension to b > and critical impact parameter 

For nonvanishing b the solutions, of course, are not axisymmetric and show a nontrivial 
azimuthal dependence on 9 = Q^ x . We shall simplify the issue by performing an azimuthal 
average on 9 and, furthermore, by performing a spin-0 projection of the a <-> a relationship. 
Instead of the simple translation x — > (6 — x) we shall take the relation (better expressed 
for the Fourier transform a(k)) 



which has no memory of the direction of b and is equivalent to an azimuthal average if a 
only depends on k 2 ,b 2 . Note that this procedure singles out the vector x corresponding 
to the source at x = with respect to the vector b — x and is therefore asymmetrical 
with respect to a and a. We can interpret it by saying that we look at the cc 2 -dependence 
of a in the average field of a, which is much similar to the collision of a point-particle 
with a ring-shaped source. Then, the symmetrical result will be obtained if we look at 
the (6 — a?) 2 -dependence of a in the average field of a. 

By replacing eq. ( 14. 1ft in eqs. (12.231) we soon realize that the ansatz a = a(x 2 = r 2 ), 
a = a{r 2 ) = (a((b — x) 2 ))g is self-consistent and, by the same manipulations as in sec. El 



(J2/E) e ~0.47, 



(3.25) 




(4.1) 
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we obtain the equations 



- (2rf)^)] = 
and, therefore, using again the function p of (13. 4ft : 

ap(r 2 ) = — -, p = 2( 7 r J R) 2 aa . (4.3) 

27T 

4.1 Solutions for 6 > and perturbative expansion 

Eqs. ( 14. 3p differ from those valid at b = by the replacement of a factor of a by a. 
In order to relate a = (a((b — x) 2 ))g to a let us note that eq. (14. ip has two distinct 
regimes, according to whether r ^> b (r <C b). In such regimes one can set, approximately, 
b = (r = 0) in the right hand side. Therefore, we are led to replace a with the simple 
approximation 

a(r 2 ) ~ a(r 2 )Q(r 2 - b 2 ) + a(b 2 )Q(b 2 - r 2 ); <j(r 2 ) ~ 6(r 2 - b 2 )a(r 2 ) . (4.4) 

Note that this approximation is exact for a® of eq. (I2.16P and for the collision with the 
ring-shaped source envisaged before. 

Introducing the above approximation in eqs. (14.3)) one has 

p(r 2 ) = ^ 6(r 2 - b 2 ) , (4.5) 

so that, for r 2 < b 2 , the "repulsive" Coulomb potential is absent and p = 0. This in turn 
leads to the solution 

p = i? 2 cosh 2 x(r 2 ), (r 2 >b 2 ); p = p(b 2 ) + p(b 2 )(r 2 -b 2 ), (r 2 <b 2 ); (4.6) 
r 2 = b 2 + R 2 (x + sinh \ c °sh X ~ Xb — sinh cosh Xb) , 

where we have introduced the hyperbolic angle x — x( r2 ; b) and the notation Xb = xip 2 )- 
The corresponding longitudinal and transverse fields are 

1 r L ' 2 dr 2 1 

<r 2 ) = - = ~(X(L 2 ) - X (r 2 )) (r > 6) (4.7) 

2tt J r 2 p[r 2 ) it 

= -(x(^ 2 )-X b ) + ^-log (n f2 2 (r<6) 
7r 27rt 6 p(0) + £ 6 H 

, i / 2\ ,ioi2j 1 _ P l-tanhx(r 2 ) n . 
= " (r > = 4191 * = W = (tt_R) 2 ' (4 ' 8) 

where £& = tanhx b , p(0) = p{b 2 ) — b 2 t b , and we have fixed the additive constant in the 
longitudinal field by requiring a(L 2 ) = 0, L being the IR cutoff parameter. 
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Because of the linear behaviour of p(r 2 ) for r 2 < b 2 , eq. (14.61) leads to the possibility 
of enforcing the boundary condition p(0) = 0, typical of the perturbative expansion, and 
required for self-consistency by the reduced action itself (see Appendix). For that to 
happen we must have 

R 2 

p(b 2 ) =R 2 cosh 2 X b = b 2 p(b 2 ) = 6 2 tanh Xb ; = t b (l - t 2 b ) , (4.9) 

a condition which resembles eq. ( 13.23!) found before. 

The criticality equation ( 14.91) is cubic in the t b parameter and determines the branches 
of possible solutions with p(0) = 0. At the critical value b 2 = 3v / 3-R 2 /2 of the impact 
parameter, the equation is stationary. For b > b c , there are two solutions with nonnegative 
p, one with t b — > 1 for b ^> b c - which will be related to the perturbative one - and the 
other with t b — > 0. The third formal solution,with t b < — 1 is actually to be discarded 
because it would require p(r 2 ) < at large distances as well. For b ^> b c , the solution with 
larger X b oo matches the perturbative solution at large distances. In fact, a simple 
iterative evaluation of eq. (14.61) yields the large-r behaviour 

R 2 4r 2 2b 2 
p(r 2 ) ~ r 2 -— log ^T^ : r 2 (xb) = R 2 ex P (l + 2x b + smh2 Xb -—^). (4.10) 
2 r\xb) R 

Since exp2x 6 ~ 4b 2 / R 2 for b ^> b c , it follows that r 2 (x fe )/4 ~ b 2 is just the scale of the 
perturbative solution in sec. \2.1\ as anticipated. 

On the other hand, for b < b c , there are no real valued solutions to eq. (14.91) with 
pip 2 ) > nor, equivalently, to the boundary condition p(0) = 0. It is not clear how to 
replace this boundary condition and thus to define a meaningful real valued solution for 
b < b c . For instance - since we cannot reach p(0) = - we can try to do our best and 
look for a Xm{b) such that p(b 2 ) — b 2 p(b 2 ) is minimal. This yields the condition 

cosh 3 x m sinhxm = t, — > ( 4 - n ) 



2R 2 " m ^ (l-t 



2 ^2 
m) 



which, for any b < b c , admits real solutions such that, while b decreases, Xm decreases from 
Xm{b c ) = Xc to Xm(0) = and p(0) increases from to R 2 . Therefore, this kind of solution 
determines Xo — as its b = limit. However, there appears to be no compelling reason 
for this choice, except perhaps that, among the real- valued solutions, the "distance" of 
this one to the p(0) = complex solution to be studied in sec. 7 is smallest. 

We remark that in this b > case, like in the problem discussed in sec. 3.2, the 
critical impact parameter b c separates - in the real-valued domain - the class of "weak- 
field" solutions having p(0) = (for b > b c ) from that of "strong-field" solutions with 
p(0) > and a small-r 0-singularity (for b < b c ). The latter solutions, however, appear 
to be somewhat ill-defined. 
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4.2 The on-shell action and its singularities 

By evaluating the action in eq. (12.221) on the solutions (14.61) we find the expression 



a(b 2 ) + a(0) - ~ J d 2 xVdVa + J d 2 x(V 2 (j)) 2 (4.12) 



2tiGs 

= 2a(b 2 ) + J a (6) + / (6) , 
where we have used a(0) = a(6 2 ) and we have evaluated the integrals 

„ 2 , 1 Z" 00 dr 2 / r 2 \ If 2v 26 2 e"» 



27r p(r 2 ) y p(^ 2 ) / 27r y i? 2 coshx 6 

In the actual evaluation, we find that the integrals (14.131) are related by an integration by 
parts yielding 1$ = —(I a + a (b 2 ))/2, except for a boundary contribution ~ p(0)(l — p(0)). 
The latter is discussed in more detail in the Appendix, where we argue that its consistent 
treatment would require p(0) = in all instances, even at the cost of picking up complex 
solutions of the field equations. Nevertheless, if we decide to keep it, we obtain, after 
simple algebra 



A(b,s) = Gs {2( X (L 2 ) - X t) + l - jp—gr- b ) ■ ( 4 - 14 ) 

It is amusing to note that the above expression for the action is stationary with respect 
to Xb at fixed b precisely when the "criticality condition" (14.91) holds. This suggests an 
alternative interpretation of the condition p(0) = 0, namely that of requiring stationarity 
in a "sum over solutions" (or perhaps better over collective coordinates contained in the 
solutions) definition of the S-matrix 

S(b,s) = J ' dxbKXb)exp(iA(b,s]Xb)) , (4-15) 

where however the integration measure p(x&)> possibly related to a functional fluctuation 
determinant around the given solution, is actually not available. 

Sticking for the moment to real- valued solutions, and using, for b > b c , the criticality 
equation in order to eliminate b 2 in terms of Xb, we finally get the convenient expression 

A(b, s) = Gs (2( X (L 2 ) - xt) + 1 - 1/h) ; (b > b c ) . (4.16) 

By evaluating Xb from ( 14.91) we find the large-6 behaviour A ~ 2Gs(log^ + Jtj), which 
checks with the perturbative expansion, the slight difference of the i? 2 /6 2 -correction being 
due to the azimuthal averaging procedure. 
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The known b- dependence of the action allows to find the elastic scattering amplitude 
A(s, q 2 ) by a Fourier transform, and the related classical deflection angle by a stationarity 
equation in b: 

^A(s,q 2 ) = ~ 2i J d 2 bexpiA(b,s) expibq (4.17) 
Qd = -V b A(b,s). (4.18) 

By applying eq. (14.181) . a simple calculation on eq. (14.141) yields the expression of the 
deflection angle 

6 c i q c i Rb R R , _ R 2 . . , , . . , _ „ . 

which determines the (resummed) corrections to the Einstein deflection as function of b, 
starting from large impact parameters down to the critical radius b = b c . Actually, the 
deflection becomes maximal (8 c i = tt) at some b > b c , the limiting value of the r.h.s. of 
(14. 19p being (4/3) > 1, showing that at such large angles the impact parameter frame- 
work is not fully meaningful. Alternatively, if we interpret A(b, s) as A(J = by/s/2,s), 
and we transform back to energy and scattering angle by convoluting A with Pj(cos6), 
we arrive at the (perhaps more physical) result: 

M& = & c ) = 2(4/3) 1/4 >7r/2, (4.20) 

meaning that, at b = b c , the two particles already invert the sign of their relative momen- 
tum. 

Note that the action (I4.16P develops a branch cut singularity at b 2 = b 2 = 6\^3 G 2 s. 
It is soon realized, using (14.161) . that, while t b has a square-root singularity at t& = t c , the 
action branch-cut is of type ~ (6 — b c ) 3 ^ 2 , because the total derivative dA/dxb vanishes 
also at Xb = Xci being sinh 2 \c = 1/2- In other words, the action is stationary in \b at 
that point, just like the criticality equation ( 14. 9ft . We expect this feature to be even more 
general than the present model, because, were the action stationary close to A c but at a 
different value, there would be another pinch of two solutions besides the one we know at 
Xb = Xc- F° r the t wo t° coincide, the action expansion around A c should start at order 
(Xb ~ Xc) 2 , the first nonanalytic piece being (xb ~ Xc) 3 - The action for b — > b+ turns out 
to have the following expansion 

= -2V3(xb-Xc) 2 + ^(Xb-Xc) 3 + 0((xb-Xcr) (4.21) 

where the analytic piece dominates around b ~ b c and provides the deflection we have 
just discussed. 
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4.3 The effective metric 



The metric describing KeS and the reduced action is obtained from eq. (I2.34p by specializ- 
ing to the axisymmetric solutions a(r 2 ) and <f>(r 2 ), which are exact at b = and averaged 
out at b > 0. We then obtain 

ds 2 = -dx + dx-(l - 2(nR) 2 Q(x + x-) ^P-) + 2ttR (a(r 2 )5(x-)(dx-) 2 + a(r 2 )5(x + )(dx + ) 2 ) 



{jRf^Lp. (\ x + \5{x~){dx~) 2 + \x-\5{x + ){dx + ) 2 ) + ds\ 



Q r 1 

ds% = (1 + 2{7cR) 2 0{x + x-)<j))dr 2 + r 2 (l + 2(<kR) 2 Q(x + x-){<P + 2r 2 '$))dQ 2 , (4.22) 

where, for b > 0, the azimuthal averaging is done at fixed x 2 = r 2 . 

We think that the effective metric so defined is really meaningful on the real valued 
lagrangian solutions for b > b c only In fact, in such a case, we have matched the solution 
with larger Xb to the perturbative expansion at large distances and, furthermore, the fields 
of eqs. (14.61) and (14. 7p are well behaved in the small-r region also. In particular, in the 
transverse part of the metric, the fields 

1 4r 2 

~ _ (0 + 2r 2 0) ~ — — log ^— (r » R) (4.23) 

o7r r r [Xb) 

have the role of decreasing the circumference over radius ratio at large distances, while 
<fi = 1 — tb becomes just a constant at small distances. 

On the other hand, if we take the real-valued solution defined above for b < b c , the 
situation does not change much at large distances but, at short distances, the fields of 
eq. (14.231) both develop a — p(0)/r 2 singularity whose interpretation is doubtful, because of 
the rescattering and string corrections neglected in the present approach. Should we take 
that singularity seriously, the coefficient of dr 2 would become negative at some r ^ R 
without any major change in the rest of the metric, because d(r 2 (p)/dr 2 is instead regular, 
except possibly at r = 0. This is perhaps one more reason to stick to the condition 
p(0) = for b < b c also, even if that means considering complex solutions (cf. sec. 7). 



5 Momentum space formulation 

5.1 Effective action and equations of motion 

In order to reformulate the problem in momentum space and to make some symmetries 
more manifest it is convenient to place the sources for the fields a and a at some generic 
points in transverse space &i and 62- Normalizing Fourier transforms as: 

a(Jfe) = ( dxe ikx a(x) a{x) = f d 2 k e - lkx a(k) = j [dk]e~ ikx a(k) , (5.1) 
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the reduced effective action (12.221) for the IR-safe polarization and after neglecting rescat- 
tering can be rewritten in momentum space as follows: 



ttA f d 2 k 



J [ e ik ^~ b ^ 7l (k) +e ik ^- b ^ l2 (k) - e ifc(b ^ b2) 7l (fc)7 2 (-fc 
nR) 2 



2 



d l k 



-h{k)h{-k) - h(-k)K(k) 



(5.2) 



Here ji(k) and h(k) are related to the Fourier transforms of a, a and h(x) by: 



7i(fe) = k ^e-^ ■ T2 (fc) = *M e -™» , h {k) = -k 2 <P{k) , (5.3) 



and 



H(k) = -^ / d 2 k x d 2 k 2 5{k -kt- fc 2 )7i(fci)72(fc 2 y (fclbl+fc2b2) sin 2 
The equations of motion that follow from this action read: 



(5.4) 



h{k) = H{k) -> fcrr(fe) for 6i -»■ 0, 6 2 -> b (5.5) 
7i (p) = l + ^p 2 j fh.h(-k)>n{k + p)e <w * sin 2 fl p/fe+p (5.6) 

72 (p) = l + ^-p 2 j ^-h{-k) l2 {k + p)e^ 2 sin 2 fl p/fe+p . (5.7) 

Eliminating h(k) through (15.51) we get two coupled equations involving just 71 and 72 : 

R 2 f d 2 k 
71 (p) = 1 + — p 2 I 7i: ^d 2 k 1 d 2 k 2 5(k + fci + fc 2 ) 



72 (p) = ! + ^ 2 



2vr 2 ^ 7 (k 2 

7i(fci)7 2 (fe 2 )7i(fe + p)^ 2 -^ sin 2 12 sin 2 # p/fc+p , 
7i(fci)7 2 (fc 2 )72(fc + p)e ifcl(bl - b2) sin 2 12 sin 2 6 p/k+p . (5. 



At least perturbatively, these equations imply the relations: 

7 2 (P) = 71(-P)=7i» , h(-p) = h*(p). (5.9) 
Setting finally b 2 = —61 = 6/2 we get the basic integral equations: 

h(-k) = Jj^ y d 2 k 1 d 2 k 2 5(k + fci + fc 2 )e i ( fc2 - fcl ) b / 2 7l (fc 1 )7*(fc 2 ) sin 2 la (5.10) 

71 (p) = l + ^-p 2 j d A- h (-k) 7l (k + p)e- ikb / 2 sm 2 e p/k+p , (5.11) 
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which can be reduced to an integral equation for 71 (p) (or 72 (p)) alone: 



R f ark 
71 (p) = 1 + — p 2 J j^d 2 kKl 2 k 2 5(k + fei + fe 2 ) 

7 1 (fc 1 )7*(fc 2 ) 7l (fc + p)e ifc2b sin 2 12 sin 2 p/fc+p . (5.12) 

It is also quite easy to check that the action, on the e.o.m. takes the simpler form: 

Aeom = ^ j^r W- lkb li{k) + e lkh l2 {k)) - 2e- fcb 7l (fc)7 2 (-fc)] . (5.13) 

An easy way to prove this is to note that our action, after elimination of <ft, is a functional 
of the 7j of the form: 

A = A X + A 2 + A 4 (5.14) 
where A n is homogeneous in 7, of degree n. By Euler's theorem: 

/X A 
dk li — = A l + 2A 2 + AA 4 . (5.15) 
hi 

On the equations of motion this combination must be zero and therefore we can eliminate 
A 4 in terms of A\ and A 2 . 

A 4 = -A 1 /4-A 2 /2, (5.16) 

which gives the claimed result. This argument only works modulo boundary (surface) 
terms. And, indeed, in the case discussed in the previous sections in which p(0) 7^ 0, such 
boundary terms can be shown to be present since their absence would imply the relation: 

It = -~[I a + a(b 2 )], (5.17) 

which only holds if p(0) = 0. There are indications that the momentum-space approach 
automatically implies such a relation. Indeed, unless there is some singularity at small 
momenta, total derivative terms are set automatically to zero by momentum conservation. 
The boundary terms due to p(0) 7^ 0, on the other hand, come from a short- distance 
boundary which does not look to be present in the momentum approach. 

In principle, the above integral equations can be solved by iteration on a computer. 
The advantage, with respect to the position-space formulation, is that the iterative solu- 
tion carries automatically with it perturbative boundary conditions. The iteration pro- 
cedure is expected to converge only for sufficiently small values of R/b where it should 
reproduce the perturbative expansion. 

Some preliminary numerical results have only been obtained [21] under the assumption 
that 7(fe) and h(k) depend only on k 2 and b 2 but not on k ■ b. This is only consistent 
with the field equations if we average them over the direction of b, i.e. if we make the 
replacement: 

e ik2b -> J (kb) , (5.18) 
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which is the momentum space version of the the azimuthal averaging procedure made in 
sect. 4. 

A first interesting indication following from the numerical analysis [21] is that, at 
sufficiently large b where the iteration converges, ^y(k) approaches unity at small k and 
a 6-dependent constant larger than unity at large k. This is fully consistent with the 
position space small-r result: 

a ^ ~ 2 * r2 , f>(b) = tanhxb = 1 - 0(R 2 /b 2 ) , (5.19) 

and confirms that, in the momentum-space approach, one has automatically incorporated 
the condition p(0) = 0. 

At some critical value of R/b the iterative solution is found not to converge any more, 
showing again the existence of a critical value for that ratio. One finds [21] (b/R) c ~ 
1.6 ± 0.1 a number that matches well (even too well!) our analytic estimate: (b/R) c ~ 
2 -i/2 3 3/4 ^ 1.6I, A more accurate numerical calculation that does not use azimuthal 
averaging appears to give [22] a slightly higher value, (b/R) c ~ 2.38. All these results are 
compatible with the CTS lower bound given in [5], i.e. {b/R)^ TS > 0.80. 

We conclude that the momentum-space approach gives numerical results for the bound- 
ary of the perturbative regime and the estimate of a critical (b/R) c that confirm those 
of the previous sections and are also compatible with the classical CTS-based collapse 
criteria. We should stress, however, that numerical momentum-space techniques, being 
based on an iterative procedure, cannot be easily extended below b = b c . Searching for 
complex solutions remains a serious challenge. 



5.2 Arguments for the existence of a critical b/R ratio 

A quick - though approximate - way to argue for the existence of a critical R/b comes 
by considering the integral equation (I5.12p for p~ 1/R < 1/b and to realize that, in this 
case, the integrals should be dominated by the regions in which all the arguments are 
roughly of the same order ~ 1/R. The equation for this "average" 7 then takes the form: 

7 = l + K 2 (R/b) V, (5.20) 

where K is a numerical constant of 0(1). It is easy to see that the perturbative solution 
7 = 1 + ... ceases to exist above a critical value of R/b, (R/b) c = 3 J^ K - Near this critical 

point the action becomes singular with a (6 — b c ) 3 ^ 2 behaviour similar to the one found in 
sect. 4. 

It is perhaps worthwhile to notice at this point an intriguing relation between eq. (15.201) 
and the equation determining the turning point r = r* for a null geodesic impinging on a 
Schwarzschild metric of radius R at impact parameter b. In this latter case the equation 
reads: 

R r* 

— = x-x 3 , x = — , (5.21) 
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giving the well-known result that the turning point disappears for b < b c = ^^-R. Simi- 
larly, eq. (I5.20p can be put in the form: 



Kj = y-y 3 , y = K{R/b)j, (5.22) 

giving, for K = 1, the same critical value for b/R. 

Amusingly, eq. ( 14. 9 p takes a similar form, this time in terms of R 2 /b 2 , i.e. 

^ = z-z 3 , z = p(b 2 ), (5.23) 

and thus gives again the critical value though for R 2 /b 2 . 

The above reasoning suggests that a more rigorous argument for the existence of a 
critical b/R could possibly be constructed along the following lines. In eq. (15.121) rescale 
all momenta by a factor |6| and distinguish the new dimensionless momenta from the old 
ones by a tilde. Also, multiply both sides of the equation by a factor R/b and define 
j3 = {R/b)^\. Then eq. (15. 12ft takes the form: 

^(fe!)/3*(fe 2 )^(fe + pie"" sin 2 9 12 sin 2 6 r/ , +r , (5.24) 

where e is the unit vector in the direction of b. 

The rhs of the new equation is a functional of j3 and a function of p. When the lhs R/b 
is very small the equation can be solved by taking (3 small so that the cubic term on the 
rhs is negligible. But when R/b is sufficiently large this is no longer the case. Taking j3 
large and of order R/b may not help to find a solution if the cubic term takes over and has 
the wrong sign. Hence it should not be impossible to show, by some functional analysis, 
that this equation does not have perturbative solutions for sufficiently large values of R/b. 
We do not attempt such a proof here. 

6 Particle production and inelastic unitarity 

So far we have neglected the imaginary part of 5, i.e. the phenomena associated with the 
production of gravitons. Even when we limit ourselves to the IR-safe (TT) polarization 
there are at least two interesting issues to be addressed: one concerns the spectrum of 
the produced gravitons (as a function of their transverse momentum); the other is the 
damping of the elastic amplitude caused by the opening of inelastic channels. We would 
like to study both effects as a function of b (or actually b/R) in order to see whether some 
interesting physics shows up as we approach a critical value. 

A convenient way to study production amplitudes within our effective action approach 
is to introduce an auxiliary source J coupled to the field h that corresponds to the physical 
TT graviton. At the same time, in order to ensure full inelastic unitarity, contributions to 
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the action from Im 5 have to be included. The standard procedure would be to do all this 
at the level of the 4-dimensional action (I2.27P by coupling the source J to the canonically 
normalized 4-D field ^V 2 $. We shall instead use a "shortcut" and modify directly the 
reduced action (I2.22p as follows: 



nA nA ttA (tiR) 2 f- 1,,,.,, IN 2zF 2VY r . . 

' d 2 k -h(k)h(-k) h h(-k) i=J{k) , (6.1) 



Gs Gs Gs 2 J y2 7T irRVGs 

where the explicit Y = log s dependence takes effectively into account longitudinal phase 
space. We will also interpret the additional terms in the action as being defined "on-shell" 
i.e. on the equations of motion of the unperturbed action. This is correct as far as the 
additional source term is concerned, but probably an oversimplification for the additional 
imaginary part in (16 .ip which presumably changes the field equations. With this caveat 
we shall now proceed to the computation of the full (inelastic) 5*-matrix. 
The functional: 

e iW(J) = e iA eom (J) ^ ( 6 _2) 

will generate, through its functional derivatives with respect to J, the scattering ampli- 
tudes for producing an arbitrary number of gravitons: 

S (2 2 + kl + k2 + . . . *„) = ( SJ( _ ki) Sn jJ( _^ ) jrf • (6.3) 

This prescription can be checked to reproduce, at lowest order, the one-TT-graviton 
production cross section -and hence the imaginary part of the elastic scattering amplitude- 
as given both by the H-diagram. In order to consider the general multi-graviton produc- 
tion amplitude we note that, at sufficiently high energies, a WKB-like approximation 
holds (since, at very high energies, S 2 W/(SJ) 2 (5W/5J) 2 , etc) : 



e iW(J) 



5J(-kt) ...5J{-k 



(^)« G5>L ■ ■ ■ (SL + subleading terms - (6 - 4) 

As a consequence, the emitted gravitons are approximately uncorrelated and the multi- 
graviton amplitude factorizes. Furthermore since, on the equations of motion, 



i8A 
5J(-k) 



= WGsYR h(k) J=0 , (6.5) 

j=o 



we find: 

5(2 _ 2 + k, + k 2 + . . . k n ) = (e tW(J) ) J=Q {iVGsYR) n J] h(h) J=0 , (6.6) 
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and 



= wniw«i-, (-) 

i 

with h(k) given in eq. ( 15. 5ft . It is easy to double check that, at the lowest level and for 
n = 1, this reproduces the one-graviton cross section discussed in sect. 2. 

At the same time, the elastic amplitude will be absorbed. Its absolute square will 
be controlled by the imaginary part of the on-shell action at J = and can be easily 
computed again in terms of h(k). We find: 

a el = \S(2 -> 2) | 2 = exp(-2/mA/=o) = exp (-GsR 2 Y J d 2 k\h(k)\ 2 ^j , (6.8) 

i.e. precisely in such a way as to ensure the inelastic unitarity of the S-matrix. 

These results can be summarized by writing the S-matrix in an operator form involving 
also the longitudinal-momentum degrees of freedom: 

S = exp (iVGsR J ^L(h(k)a k + h(k)* = ex P (-GsR 2 J 7^IM*0| 2 ) 

exp (iVGsR J ^h(k)*a£j exp (iy/GsR J ^Lh(k)a k ^J , (6.9) 

where , a\ are canonically-normalized destruction and creation operators of physical 
gravitons of momentum k and IR-safe polarization and we have left out for simplicity an 
overall c-number phase exp iA(b, s) containing the on-shell uncorrected action. This gives 
back, for instance, eq. (16.81) after realizing that the longitudinal momentum integration 
just provides a factor Y = logs. 

The S"-matrix (16.91) when acting on the Fock vacuum of the a& , a\ operators creates 
a coherent state of physical gravitons in which we can compute the expectation value of 
the associated canonical quantum field that we denote by h can . This is best done by using 
the LSZ formalism according to which: 

^-^e lkx d^(2\h can (x)\2)\^ , (6.10) 
where, in our case, the canonical TT-graviton field is given by: 



> tn 



/>can = (SttG)- 1 / 2 -^ , (6.11) 

and the extra factor l/y/2 comes from our normalization of the polarization tensors. 

Using properties of the coherent state generated by ( 16.91) it is easy to check that the 
following metric fluctuation satisfies ( 16. lOf) : 
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where, as previously explained, eJJ act as differential operators and we have inserted an 
ie prescription although the LSZ formula is only sensitive to the principal part of the 
propagator. 

Equation (16.121) should be consistent with the effective metric of eq. (I2.34p and, at 
lowest order in R/b, with the one that follows from eq. (12.141) . Indeed, if one looks at the 



principal part contribution, one finds, in the small- fc limit 1 implicit in our procedure: 

(h(x), w ) = (ttR) 2 J d[k]e T Jh(k)&(x + x~)e~ lkx = (7rR) 2 e T J&(x + x-)V 2 <f) , (6.13) 

which explains the normalization of the <fi field used in eq. (12.341) . 

Using the small- A; limit of h(k) we find h c i(x) ~ R 2 /r 2 at large r, in agreement with 
(I2.14p but in apparent disagreement with the standard quadrupole formula, which would 
require h c[ (x) ~ R 2 /br. This can possibly be explained by the fact that the TT polariza- 
tion is emitted mainly in the forward and backward direction \k\ 3> \k\ so that there is 
still a non-trivial flux of TT-gravitons at null infinity. On the other hand, the (IR-unsafe) 
polarization that we neglected exibits a 1/r behaviour in agreement with the fact that IR 
singularities are associated with classical radiation. 

It is very tempting, at this point, to guess a generalization of our result (16. 9p to include 
the IR-sensitive polarization. This would read: 

S = exp ^\/GsR J (ih TT (k)a k — h LT {k)b k — h.c^j , (6.14) 

where Htt — h and: 

h LT = -^ [ d 2 k l( i 2 k 2 5{k + kx + fc 2 )e i ( fc2 - fcl ) b / 2 7 1 (fci)72(fc 2 ) sinfl 12 cos# 12 . (6.15) 
n z k z J 

This would provide a unitary S-matrix whose matrix elements, however, are only finite 
within particular coherent states that include soft bremmstrahlung, a well-know situation 
in perturbative QED and quantum gravity, already discussed in [3] . Once this is properly 
done there should be no major conceptual obstacle in including the effects of the LT 
polarization on the gravitational collapse problem. 

What remains to be done is to evaluate h(k) in different regimes in order to extract 
both the spectrum of the emitted gravitons and the absorption of the elastic amplitude as 
a function of b/R. This can only be done, of course, after solving the classical equations, 
similarly to what already done for Re^4 in position space. We can summarize our present 
understanding on this matter as follows: 



5 Keeping the fc-dependence amounts to multiplying the O-function by Jo(\k\\/ x + x ), which implies 
the cutoff x + x~ < R 2 mentioned in sec. 2. 3. 
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• At b ^> R the spectra have a logarithmic behaviour at small k and an exponential 
damping at k ^> b^ 1 . In other words the typical tranverse momenta of the produced 
gravitons are, not surprisingly, of the same order as those of the exchanged gravitons 
(see the discussion in sec. 2). In turn, the elastic amplitude is suppressed by 
exp (— GsYR 2 /b 2 ). Note that, for R > l s , such a suppression dominates over the 
one due to string excitation, that we have neglected. 

• As one approaches the region b — > b c ~ R from above, physics appears to be 
rather smooth. Nevertheless, the graviton spectrum is now cutoff at momenta of 
order 1/R (i.e. of order of the Hawking temperature of a BH of mass and the 
elastic amplitude is suppressed by an exponential factor exp (—cGsY) (with c some 
constant of 0(1)) which, modulo the factor Y = logs corresponds to exp (S), with 
S the Bekenstein-Hawking BH entropy 

• It would be very interesting to find out what happens if one goes to the region 
b << R, in particular whether the cut-off on momenta keeps growing like 1/6 or 
remains "frozen" at 1/R as black-hole evaporation would suggest. It is not yet fully 
clear how that region can be studied. The possibility of an analytic continuation of 
S-matrix elements and effective action solutions for b < b c is discussed in the next 
Section. 

7 Action and complex solutions for b < b c 

We have realized in sec. H] that real- valued field solutions with p(0) = exist only for 
b > b c , and that, below b c , they become complex. We have also shown that the action has 
a branch-point singularity at b = b c , presumably due to the pinch of two such solutions 
of the criticality equation (14. 9ft . The problem then arises of how to define both S-matrix 
and solutions for b <b c . 

It is tempting to try the simplest possibility, and to continue analytically the S-matrix 
on the basis of the form ~ (b 2 — 6v^3G 2 s) 3 / 2 of the branch-cut, by choosing a physical 
energy-sheet reached by an s + ie prescription. According to the expansion (14.211) this one 
corresponds to the Xb~Xc determination having a negative imaginary part and contributes 
a positive imaginary part to the action. This criterion tells us to take that particular 
complex solution as the physical one on which the action, and thus the S-matrix, should 
be computed. 

An additional argument for this choice comes from the tentative interpretation of the 
p(0) = condition as stationarity equation of the integral over Xb in eq. (14. 15ft . If we 
define the analytic continuation of the action by that integral - which, for b < b c has 
two complex conjugate stationarity points - we should take the one for which the saddle 
point is stable. Then, by ignoring the measure factor, calling Xb the stationarity point 
and simply \ the integration variable, we have the expansion 



•A(b,s;x) -A(b,s;xb) 
Gs 



2(-l + 



b 2 sinh Xb 



)(X-Xb) 



(7.1) 



R 2 cosh 3 Xb 
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+ — (1 - t 2 b )(l - 3tl)(x - Xb f + 0(( X ~ Xbf) ■ 

We note again that the action is stationary on the solutions of (14.91) and that the fluctu- 
ation coefficient in eq. (17. f p becomes, for small b — b c , 



i(l - 3t 2 b ) ~ T2V2t^j 1 - (t b - to) * ±i 1 - p • ( 7 - 2 ) 

The saddle point is then stable when the above coefficient is positive imaginary, corre- 
sponding to damped fluctuations, yielding again the solution with Imt^ < for which the 
action acquires a positive imaginary part, as noticed before. 

By then taking for b < b c the complex solution for tf, or Xb with negative imaginary 
part, we obtain, from eq. (14. 2f p . 



IH^M = ^ (1 _| P , |S(M ^ exp (_^ (1 _| Pl . ,, 3) 

Equation (17. 3p implies that the S-matrix, when analytically continued to b < b c , acquires 
an additional absorptive part (on top of the one due to TT-graviton production discussed 
in the previous section) whose interpretation calls for the opening up of some extra chan- 
nels in this new regime. It is tempting to think of these as quantum analogs of the 
black- holes that are expected to be formed on glassical grounds [5J. At the same time, 
the dominant solution is also complex-valued, with Imp(r 2 ) < for r ~ R. 

The above features are confirmed by evolving the complex solution to smaller values 
of b. In the b — > limit the stable determination becomes, by eqs. (14. 9p and (17.21) . 

t b ^e-^y\ X, = -^ + e- /3 (^) 2/3 , (b«R), (7.4) 

while the complex conjugate solution is unstable. Correspondingly, the action for b <^ R 
takes the form 

A{b,s) = G S (2 X (L 2 )+m-3e^ 3 (-|) 2/3 ), \S el (b,s\ 2 ~e- R ^-^W 2/3 K (7.5) 

R 



For 6 = the suppression factor in the elastic cross section is just exp(— iSbh(\A/2))? 
with <Sbh(M) the Bekenstein-Hawking entropy of a Schwarzschild black hole of mass M. 
This is compatible with a statistical interpretation where a fraction 1 / a/2 of the incoming 
energy goes into forming a black hole, even if the relationship of our S-matrix framework 
to such a semi-classical statistical picture is yet to be clarified. 

On the other hand, the suppression (17. 3p of the elastic amplitude appears to die off as 
(I — p) 3//2 for b — » b~ . This looks like an interesting (and we believe robust) result calling 
for a physical interpretation. From a classical standpoint, this limit should correspond to 
the production of a nearly extremal Kerr black hole with J < J c = GM 2 . If the mass 
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of the produced black hole would remain finite in this limit its entropy would approach 
a finite value (just half of that of a Schwarzschild black hole of the same mass) and the 
statistical interpretation invoked earlier for b ~ would fail. However, it is conceivable 
(although, to the best of our knowledge, still not proven) that in a nearly critical-collapse 
situation most of the initial energy and angular momentum are radiated away to infinity 
leaving only a vanishing mass and angular momentum to collapse at the critical point. 
In this case, if we insist on identifying our elastic suppression with an entropy factor, we 
have to assume that the mass of the Kerr black hole being formed in the collision vanishes 
like (b c — 6) 3 / 4 for b — > b~ . This would imply a "Choptuik exponent" of 0.75 in our critical 
collapse, i.e. about twice the original exponent of ~ 0.37 found in Choptuik's original 
paper g6j@ 

Corresponding to the small-6 parameters (17.41) . the expression for the physical field 
solution p(r 2 ) has the initial value 

^1 = —t h = e ^(-) 4/3 (7 6) 

which is consistent with p(0) = and has a small, positive real part also. This is sufficient 
to have an r 2 -evolution of p(r 2 ) with increasing real part and negative imaginary part, 
which tunnels to a perturbative, real valued behaviour at large distances. We can further 
check that the action integrals in (I4.13P are well defined on the physical solution, and their 
evaluation holds unchanged, except that p(0) = is now built in, so that the result (14.141) . 
and thus (I7.3P and (17.51) obtain automatically, with the appropriate (complex) values of 
Xb and t b . 

In order to better understand the r 2 -evolution, it is convenient to come back to the 
axisymmetric b = case, by looking for possibly complex solutions satisfying the boundary 
condition p(0) = 0. Since the parametric expressions in eqs. ( 13.91) and ( 13.101) are still valid, 
we obtain p(0) = by setting xo — —in/2. With this boundary value, for r <C R it is 
convenient to look for solutions x(r 2 ) = —mj2 + r](r 2 ) such that, by ( 13.91) . 

r 2 4 
2— = 2x + sinh 2x + iir = 2r) - sinh 2r] ~ — rf + 0{rf) . (7.7) 
R 2 3 

This equation admits in turn three branches, according to the values e = (e ±l7r//3 , —1) of 
the three cubic roots of (—1), as follows 

/ 9n / 3r 2 . i p(r 2 ) 9 3r 2 2 . , 

V(r 2 ) * e(^)a, ty. ~ - v 2 ~ e*( w )s, (r<J2) (7.8) 



R 2 * - -WZ^r 2/3 , (2-r) 2 ^ - £ 



2tt y 2R 2J ' v ' r R 2 v 2R 2 ' 

We thus see that the solution with e = e l7T ^ 3 matches the physical one discussed above in 
the b — > limit. We also see that the small-r 2 exponents are fractional, and the delta- 
function flux from r = in (13.31) is shared between a and <fi. Note that there is also a 



6 One of us (G.V.) would like to thank Steve Giddings and Don Marolf for interesting discussions about 
this possibility. 
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real- valued solution (the one with p < noticed before), which however is unphysical: it is 
quite nonperturbative (p = — 1 in the large-r limit) and yields a quadratic IR divergence 
in the action. Note finally that the large-r 2 behaviour of the physical fields is instead 
perturbative, the value of xo = — wr/2 contributing (via the scale f 2 in (13.111) ) some 
subleading imaginary part. 

For b > 0, we note that the small-r behaviour of <ft changes in the r < b region. In 
fact, since in this region p = 0, then p = % must be a constant and, due to p(0) = 0, <fi 
will be a constant too: 

(27ri?) 2 = 1 - i b = 1 - e* (^)~ 1/3 ; = (r<6). (7.9) 

Therefore, b > acts effectively as a cutoff for the small-distance behavior cf) ~ d(r 2 <fi) / dr 2 ~ 
r -2 / 3 occurring at 6 = 0. The field becomes completely regular. 

We thus see that the small-r behaviour of the complex physical solution has now 
changed - compared to the real-valued ones - to a "weak-field" profil^j], due essentially 
to the p(0) = property. The latter condition acts in this context as a quantization rule, 
yielding a well-defined solution whose classical counterpart, if any, one should classify as 
being "untrapped". Therefore, at quantum level, the physical solutions for b < b c show 
no evidence of a field being confined behind the would-be horizon. 

The main challenge remains, of course, the interpretation of such complex solutions 
and of the extra absorption found in the elastic channel. A possible way to proceed is to 
continue below b c the production amplitudes of the multi-particle channels considered so 
far, and look for unitarity integrals that could make up for the extra absorption. Since 
the field h ~ 1 — p acquires an imaginary part, the latter could perhaps be interpreted 
according to the /i-definition in sec. (2.1): the imaginary part would then simply be 
the LT polarization which, though excluded in the beginning, is turned on necessarily 
below b c . In any case, within this interpretation, the Hawking evaporation required by 
the additional absorptive part should be looked for in the various contributions of such 
imaginary part to the multi-graviton spectra. Since there is no ~ log s longitudinal phase 
space enhancement, it should be quite central, with fc 3 ~ |fe| ~ 1/R and an emission yield 
~ Im^4. 

Another possible interpretation (see sect. 7) is that the missing probability goes into 
the formation of some new bound states. One may object that such bound states, be- 
ing very massive, should decay into light particles (gravitons in our case), behave as 
resonances, and that consequently the 5-matrix should already be unitary in the multi- 
graviton Hilbert space. However, it is quite possible that, within our approximation that 
neglects corrections of relative order h/Gs = Mp/s <C 1, such resonances are actually 
stable and have to be included in the possible final states in order to ensure unitarity. 
For instance, according to standard lore, black-holes have a lifetime of order R and thus 
propagators of the type (s — M\ H +iMp)~ l . If objects of this kind, even if not necessarily 

7 We use inverted commas here since the perturbation of the metric is actually still large at r ~ R and 
even more so at r ~ b. 
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to be identified as black holes, were responsible for the extra absorption below b c , their 
finite-width effects would be lost in our semiclassical regime, Gs/h 3> 1. 

8 Summary and outlook 

Let us summarize the method we have used, and the main assumptions and results of our 
investigation. 

Working in the superstring approach to scattering amplitudes in the transplanckian 
regime Gs S> H [HOE], we have used the effective action framework justified in |4] for the 
case R, b 3> A s , in which string-size effects are not very important. We have then neglected 
the so-called rescattering terms (as it should be justified for fixed (not too small) b/R) 
and, in order to avoid the known |3J but technical treatment of the IR problem, we have 
considered the emission of only one graviton polarization, the IR-safe one, described here 
by the scalar field 0. Both are technical simplifications that could in principle be waived 
one by one. 

In the framework just described (sec. 2), the dependence of the effective action and 
fields on the longitudinal coordinates simplifies, so that the problem reduces to a trans- 
verse two-dimensional effective action, that we have considered in both configuration 
(sees. 3, 4) and momentum space (sees. 5, 6). In momentum space, an iterative procedure 
for solving the equations of motion has been set up, and is suitable for numerical com- 
putations. In configuration space, the equations are studied in the axisymmetric case, in 
which they reduce to ordinary differential equations. Of course, axisymmetric solutions 
are directly relevant only in the 6 = limit, while for b > they imply the azimuthal 
averaging procedure explained in sec. 4. 

A key point we have discussed throughout the paper is about the boundary conditions 
to be set in order to determine both field solutions and action. One is provided by 
matching the perturbative behaviour for r 3> b,R. The other (p(0) = 0) is expressed 
in terms of the auxiliary field p = r 2 (l — (27iR) 2 <j)) and is still suggested by the weak 
coupling regime valid for b ^> R. In addition, we argue that it has to be valid in the non- 
perturbative regime also, as the only consistent way to treat the r = boundary. Under 
such conditions, we are able to provide analytic field solutions, the corresponding action 
and effective metric and then the phaseshift operator resumming the R 2 /b 2 corrections to 
the eikonal and the S-matrix. 

The perturbative resummation diverges at a critical value of the impact parameter 
b = b c ~ R = 2G-\/s, which separates the class of real-valued (b > b c , sec. 4) and complex- 
valued (b < b c , sec. 7) solutions, all satisfying the boundary condition p(0) = which plays 
the role of quantization condition of the problem. We also find that our estimates of the 
b c j R ratio are compatible with the classical lower bound for CTS formation, suggesting 
that our non-perturbative regime is likely to be in correspondence to classical collapse. 

For b > b c , the S-matrix has essentially the form of a unitary coherent state operator 
from which elastic absorption and inelastic spectra can be computed. For b ~ b c , the 
graviton spectrum is cutoff at transverse momenta of order h/R = Th, the Hawking 
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temperature of a black hole of mass y/s, and the corresponding elastic absorption is 
~ exp(— const. GsY). 

For b < b c , the analytically continued physical field solution with p(0) = is com- 
plex, and yields an additional elastic absorption compared to that just mentioned. The 
absorptive suppression is exponential, and the exponent vanishes like Gs(l —b 2 /b 2 c ) 3 ^ 2 for 
b — > b c and, for b = 0, is just 2-nGs. This behaviour is compatible, as order of magnitude, 
with the entropy of a black hole of properly chosen mass, as argued in sec. 7, even if the 
relationship of our S-matrix coherent state to such a classical object is yet to be clarified. 

Our field solutions provide an effective metric also, which is of shock wave type for 
the longitudinal fields, and of finite wavefront for the mostly transverse one. For b > b c , 
the profile function of the transverse field is everywhere regular in the transverse r 
coordinate. For b <b c some metric components become complex and their interpretation 
is open to discussion but, surprisingly, the <fi field keeps being regular for < b < b c also, 
and has a mild fractional r = singularity for b = 0. We refer to this feature as a "weak 
field" situation (see however our footnote in sect. 7). In this sense, our physical solutions, 
for any value of b, show no evidence of a field being confined in the small-r region. 

For b < b c , we have studied real valued field solutions also, which exist at the expense 
of violating the boundary condition p(0) = 0. They have in fact a positive p(0) and 
thus show a ~ — p(0)/r 2 singularity and a "strong field" situation. Such solutions are 
however ill-defined (depending on the value of an arbitrary parameter), just because the 
p(0) = condition is not met. For this reason, we believe them to be unphysical. Then, if 
the above "quantization" condition is imposed, the solutions for b < b c become complex 
and change to the "weak field" profile discussed before. The correct solution is therefore 
no longer confined, in a way suggestive of a quantum tunnel effect. 

The additional absorption found before for b < b c calls for extra production channels 
for the S-matrix to be unitary, on whose nature we have made a couple of guesses in sec. 7. 
Here, further work is needed in order to continue the appropriate production amplitudes 
below b c and thus to check whether inelastic unitarity is really verified: in this respect, 
our results are only partial. Nevertheless, we feel that some new physics is emerging in 
this simplified, but consistent quantum-gravity treatment that we have proposed. The 
picture outlined by our results suggests that this is the right framework for at least asking 
the questions, even if we only have some of the answers. 
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Appendix: Boundary terms and the p(Q) = condition 

Here we investigate the boundary terms present in the definition (I2.22p with the purpose 
of understanding in a more formal way the boundary condition p(0) = 0. The expression 
( I2.22p is nominally of fourth order in the derivatives of 0, but can be reduced to second 
order by introducing as fundamental field a first derivative. Furthermore, the expression 
depends on the current 7i which, for generic fields a and a, is nonlocal. We shall thus 
look for a local form of the effective lagrangian, by specializing for simplicity to the 
axisymmetric case, in which all fields are functions of r 2 only. 

We thus express the <fi field in eq. (12.221) in terms of the p function introduced in 
eq. fl4.3p .namely 

p(r 2 ) = r 2 (l-(2vr J R) 2 0(r 2 )), (A.l) 

and we solve for T~C(r 2 ) in terms of ad by the same manipulations of sees. 3 and 4. We 
then obtain 

A ~ a(b 2 ) + a(0) + 2tt f°° dr 2 ( -r 2 dd + 2H(r 2 )(l - p) - ^"^H ; (A.2) 



2nGs v ' v ' J Q V. (2ttRY 

POD 

2H{r 2 ) = / dr 2 d(r 2 )d(r 2 ) , (A.3) 



where we have specialized to the determination of 7i which vanishes at large distances - 
similarly to what we have done for a on the equations of motion - in order to recover the 
perturbative behaviour in that region. We then perform an integration by parts in order 
to eliminate the r 2 -integral in the expression of 7i 

POD 

dr 2 27i(r 2 )(l - p) - d(r 2 )a(r 2 )(r 2 - p(r 2 )) = p(0) / dr 2 d(r 2 )a(r 2 ) , (A.4) 



and we thus obtain, in the right hand side, the boundary term we were looking for. Note 
that the latter is strictly speaking nonlocal as well, because it couples p(0) to all values 
of a(r 2 ). So, if we require locality, we must set p(0) = 

Nevertheless, if we decide to keep the boundary term and the p(0) parameter, by 
replacing (1A.4I) in fl A. 2|) we get a more conventional form of the action 

" 4 - a(b 2 ) + a(0) + 2tt [°° dr 2 ( -d(r 2 )d(r 2 ) (p(r 2 ) - p(0)) - ^~^ ) . (A.5) 



2txGs y ' w J V ( 27rR ) 

We note at this point that the form (IA.5I) of the action yields the same lagrangian equa- 
tions as in sees. 3 and 4 but with p replaced by p = p{r 2 ) — p(0), which must vanish at 
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the origin. This is yet another reason for this boundary condition, which can be met by 
real solutions for b > b c and by complex ones for b < b c . 

Furthermore, the action is functional of p only and, evaluated on the equations of 
motion takes the form 

A(b, s) = 2nGs(a(b 2 ) + ^±^!2 ) = 2nGs(a(b 2 ) - 1+) = Gs(2( X (L 2 ) - Xb ) + 1 - i) , 

l tb 

(A.6) 

in agreement with eq. ( 14. 16ft and with the momentum space relationships. Thus, in this 
alternative point of view, p(0) is an additive constant in the definition of the field r 2 4> in 
terms of h ~ (1 — p), which appears in the effective metric but not in the action. Since 
this constant is not there in the perturbative regime b > R, it should finally be absent 
altogether. 
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